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A GENERAL FORM OF THE CO- MOVING TENSORIAL DERIVATIVE 


By Earle K. Huckins III and Richard E. Turner 
Langley Research Center 

SUMMARY 

A general expression for the co- moving derivative of a tensor is derived. A vari- 
able describing the coordinate velocity field is introduced. Time dependency of the 
metric elements is expressed in terms of this velocity field. The resulting description 
of motion is one of which the Eulerian and Lagrangian viewpoints are special cases. This 
general description is useful in problems involving moving boundaries or discontinuities. 

INTRODUCTION 

The field equations (as opposed to the constitutive equations) of any continuum are 
simply statements of the conservation of mass, momentum, and energy. These equations 
can be written in a tensorial form (ref. 1) which is, by design, valid in any classical coor- 
dinate system. The most common practice is to adopt an inertially fixed coordinate sys- 
tem (Eulerian formulation) in which to describe the motion. Another classical, but sel- 
dom used, approach is to fix coordinate points in identified particles of the continuum 
(Lagrangian formulation, ref. 2). In most problems, the Eulerian approach is preferable 
and completely satisfactory. However, in problems which contain moving boundaries, 
free surfaces, or discontinuities, finite difference algorithms developed from an Eulerian 
description become extremely complex at these boundaries. On the other hand, the 
Lagrangian approach, which overcomes the problems associated with moving boundaries, 
can lead to extremely complex grid networks in the presence of vorticity (ref. 3). Both of 
these problems suggest the desirability of a mixed Eulerian- Lagrangian description which 
can be obtained by introducing a completely general coordinate system. Time dependence 
of the space metric introduces additional terms in the governing differential equations. 

The equations of motion for surface flow are derived in reference 1 by considering time 
dependent metric elements. By a similar approach, the methods can be extended to three- 
dimensional problems. 

The purpose of the present paper is to derive an equation for the tensorial co-moving 
derivative of a tensor. In order to define the arbitrary coordinate system, a variable 
equal to the velocity of coordinate points relative to the primary inertial system and 
expressed in the arbitrary coordinate system, is introduced. Differential equations in 
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time for the elements of the metric are expressed in terms of this variable. By using 
this approach, it is possible to develop the field equations in any arbitrary classical 
coordinate system. 
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SYMBOLS 

tensor quantity 

components of a general tensor A of any order 

denotes that tensor A is resolved in y s 
mean divergence of coordinate velocity field 
i,j element of the metric of y s 
identity matrix 

denotes evaluation in inertial coordinate system 
surface normal 
scalar invariant 

velocity of the nodal points of y s in primary inertial coordinate system 
denotes covariant derivative of S 1 with respect to yj 
matrix formed by elements T\ j 

coefficient of covariant transformation law (see eqs. (7) and (8)) 
time 

velocity of continuum relative to x s 
control volume for defining mean divergence 
velocity of continuum relative to y s 
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X s primary inertial rectangular Cartesian coordinate system 

y s arbitrary coordinate system 

Sj Kronecker delta 

e.j k permutation symbol 

X inverse of matrix T, X = T~* 

X^ i,j element of the matrix X 

<f* coordinate velocity potential 

fi 1 angular velocity tensor 

Repeated indices denote summation over the range of the index; Roman indices 
(for example, i,j) denote a range of three; and Greek indices (for example, a,0) denote 
a range of four. 


ANALYSIS 

In the analysis to follow, the co-moving derivative of a tensor quantity is defined, 
and physically described, in an inertial rectangular Cartesian reference frame. Based on 
this definition, an expression for the co-moving derivative in a completely arbitrary coor- 
dinate system is formulated. 

Description of Co-Moving Differentiation 

Consider the motion and deformation of a continuum in an inertial rectangular 
Cartesian coordinate system x s . Let A^*'™ be the components of a tensor quantity A 
representing a characteristic of the continuum. Assume that the components A^’"™ 
have a functional form given by 
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.n 


(x s , t) 


( 1 ) 


The co-moving derivative of a variable is defined (ref. 2) as the total time derivative of 
that variable evaluated for an identified particle in the continuum. That is, the co-moving 
derivative, as the name implies, is the time derivative seen by an observer moving with 
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the particle and is evaluated in an inertial frame. Assuming that ‘ ‘ ‘ m has continu- 

Py. . .n 

ous first partial derivatives, the total time derivative can be expanded by using the chain 
rule of differentiation and written 


dA 1J - m aA 1] * ,,m 
P q--n - pq--n , 

dt at 


rwij...m 

8 pq...n dx s 


a x u 


dt 


which is equivalent to 
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_P^R + 

at 


dA ib..m 

pq- n 

9x s 


( 2 ) 


( 3 ) 


where u s are the rectangular Cartesian components of the continuum velocity in the 
primary inertial frame. The first and second terms on the right-hand side of equation (3) 
are commonly called the local and convective derivatives, respectively. This form of the 
co-moving derivative describes the Eulerian viewpoint of change in A and is of funda- 
mental importance in physical and mechanical laws. 

General Form of the Co-Moving Derivative 

In many problems, it may be advantageous (as will be discussed later) to express 
physical laws (and hence, the co-moving derivative) in a completely arbitrary coordinate 
system. In such instances, the form of the co-moving derivative must be expressed in 
tensorial form. In the following development, a general form of the co-moving derivative 
is derived for a covariant first-order tensor. These results are then generalized to mixed 
tensors of any order. 

Consider a completely arbitrary coordinate system y s superimposed on the inertial 
space. Typical parametric line segments for such a coordinate system are shown in fig- 
ure 1. Consider also an identified particle of the continuum moving about in this space. 
Assume that associated with this particle is a first-order covariant tensor A whose 
components A s have continuous first partial derivatives. Let A s (x) and A s (y) 
denote covariant components of A in the x s and y s coordinate systems, respectively. 
Treating time as an invariant, the complete transformation (including time) between the 
components can be written 

V y)= ^5V x) < 4 > 
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ARBITRARY COORDINATE SYSTEM 



where 



This transformation can also be written in the partitioned form 


(5) 



where 


T- • 
i] 



( 6 ) 


(7) 


and Roman indices have a range of three and denote spatial components. Therefore, 


Aj(y) = T.j A.(x) 


( 8 ) 


That is, the spatial components (covariant) of A transform independently of the fourth 
component. Note, however, that the fourth component does not transform by invariance. 
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Now consider the contravariant components of A. These components transform 
according to 



which can be written 




where | y s indicates y g was held constant during the partial differentation. Note also 

from equation (10) that the fourth covariant component transforms by invariance. However, 
the spatial components transform independently of the fourth component if, and only if, 
either the fourth component is identically zero or the coordinate transformation is time 
independent. Since time dependence of the transformation is an essential element in the 
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present analysis and a restriction to spatial components is desirable, the analysis will be 
limited to those tensor quantities whose fourth component is identically zero. This 
restriction eliminates consideration of velocities relative to moving coordinate systems. 
However, velocity relative to the primary inertial frame and force vectors are admissable 
and of primary importance in formulation of the field equations. 

Based on the above discussion, the transformation of the spatial covariant compo- 
nents of A can be written 


Aj(x) = A.(y) 


(14) 


Since (for a time invariant coordinate transformation) time differentiation does not alter 
tensor character, the co-moving derivative of A^x) is a tensor and can be expanded in 
the form 


dA i (x) 

dt 




9 A. (y) 

J i ir 

8A (,) 

+ A jW 

!^i + v k !iii' 

St 

>> 

CD 

St 3y k 


( 15 ) 


where the 


v k are defined by v k 


dy k 

dt 


and are the contravariant components of the 


continuum velocity relative to the arbitrary coordinate system. The components of this 
derivative in the arbitrary coordinate system are given by 


[dA.(y)' 


dt 


T u 


dAj(x) 

dt 


(16) 


where ( )j indicates that the time differentiation was performed in an inertial coordinate 
system. Substituting equation (15) into equation (16) gives 



8A i (y) 9A.(y) 


at 


+ V 1 


8y“ 


+ T U 




k o 8X ik 
_JS + vS^LS 

at 8 y s 


A t (y) 


(17) 
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Equation (17) can be rearranged and written as 



0 A.(y) 


at 


+ V" 


0 A.(y) 

0 y s 


9\.. 

+ T. . — ^ A k (y) 
ay 


3X„ 


+ T- . — A, (y) 
l] at k v,/ 


(18) 


As is shown in appendix A, the bracketed term in equation (18) is the covariant derivative 
of the quantity A^(y). Therefore, 


dAj(y)\ 

r~dt /j 


9A i (y) 


at 


+ v k 


A. (y) . + T- • 
i u; J s T 


ax. 


jk 


at 


A t (y) 


(19) 


The last term in equation (19) accounts for time dependence in the metric of the 
arbitrary coordinate system. Since, by definition 


T u x ik = 1 


then, 


T,, 


8X„. 9T„ 


ij at at jk 


Also, since 


T. • = 
i] 


9x1 

9 y* 


then, 


9T. . 


at 



( 20 ) 


( 21 ) 


( 22 ) 


(23) 
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Therefore, 


T.. 

i] 




^ 9 x-* 

8 1 

a y i‘ 

v 8t 


X, 


jk 


(24) 


Since the arbitrary coordinate system must be specified in some manner, a new 
variable S is introduced and is defined as the velocity of coordinate points of the arbi- 
trary coordinate system relative to the primary inertial coordinate system. In the arbi- 
trary coordinate system, the contravariant components of S are 


S 1 = - 


V 

at 


(25) 


Therefore, 


8x J 

at 


a y s at 


= T • S' 
S] 


(26) 


and substituting equation (26) into equation (24) results in: 


0X. 


i] at 


9S k 9 T . 

s s x 

3y* 


jk 


(27) 


Again, as is shown in appendix A, the term on the right-hand side of equation (27) is the 
covariant derivative of the coordinate velocity; that is, 


T-. 

i] 




(28) 


Therefore, by substituting equation (28) into equation (19), the co-moving derivative 
becomes 
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( 29 ) 


fdA.(y)\ 

9A.(y) . 


u i- 

+ v s A. (y) 

at i w | B 

- S k (i A k (y) 


In order to carry out the covariant differentiation in equation (29), elements of the 
metric must be determined. Elements of the metric can be related to the coordinate 
velocities through a differential equation as follows. The required metric elements of 
the y coordinate system can be written 



Sy 1 3yi 


T ik T jk 


Differentiating gives 






From equations (21) and (28), 



which can be written 




Substituting equation (33) into equation (31) and reducing gives 


(30) 


(31) 


(32) 


(33) 


3g. . 
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9t 


s ip sP |j + 8 


S q |. 
jq li 


which is equivalent to 


ag ii 

— ^ - S.i. 

0t *|j 



(34) 


(35) 
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The equivalent form of the co-moving derivative (eq. (29)) for a first-order con- 
travariant tensor can be formed as follows. Let Q be a scalar invariant. For a scalar, 
which transforms by invariance, the co-moving derivative (eq. (3)) reduces to 


dQ 

dt 



Consider a particular scalar generated by the inner product 


(36) 


Q = A l (y) Aj(y) (37) 

By generating the co-moving derivative of Q by product differentiation and substituting 
from equation (29) for the time derivative of the covariant components of A, the co-moving 
derivative for the contravariant form is found to be 



The expressions for the co-moving derivatives (eqs. (29) and (38)) can be generalized 
to any order tensor by a process of inner product formation and product differentiation. 

For instance, let 



(39) 


(40) 


The time derivatives of A J and B. can be expanded by using equations (29) and (38). 
From this result, the co-moving derivative of a second-order covariant tensor is found to 
be 
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By using a similar procedure, the results can be completely generalized to a mixed 
tensor of any order. This result can be written 


( 42 ) 


RESULTS AND DISCUSSION 

In the discussion to follow, the utility of the generalized co- moving derivative is 
demonstrated. Eulerian and Lagrangian descriptions are shown to be specialized cases 
of the generalized form. The significance of coordinate velocity derivatives for a rotating 
Cartesian coordinate system is discussed. Finally, application of the generalized form is 
described. 

The Eulerian viewpoint is generated by requiring the arbitrary coordinate system to 
be time invariant. That is, let 

S 1 = 0 ( 43 ) 

Under these assumptions the time derivative is evaluated at an inertially fixed point and 
the co- moving derivative reduces to 





(44) 


In the Lagrangian formulation, coordinate points in the arbitrary coordinate system 
follow particles of the continuum. This viewpoint can be generated by assuming that the 
continuum is stationary in the arbitrary coordinate system 
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(45) 


V s 5 0 

and that the coordinate points move at the continuum velocity u 1 ; that is, 


S* 


= u* 


(46) 


Note that the continuum velocity (and therefore, the coordinate velocity) is resolved in the 
arbitrary coordinate system. Under these assumptions, the co-moving derivative 
becomes 


+ -V ♦ “V <•::» ♦ • 


. n m i (v) 

r A pq. ..n' y/ 


u |p A sq...n W U | q A ps...n W 




(41) 


As an additional example, let the arbitrary coordinate system be a rotating rectan- 
gular Cartesian coordinate system. In order to examine the significance of the coordinate 
velocity, assume a rigid continuum with body fixed axes such that 

V s = 0 (48) 

Consider the co-moving derivative of the velocity duVdt at a point (y 1 , y2, y3) in the 
body axes system. For this situation the co-moving derivative becomes 


du^ 

dt 




where 



(49) 


(50) 


Since the arbitrary coordinate system is Cartesian, 



(51) 
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Therefore, the co-moving derivative becomes 


du 1 _ 8U 1 | as 1 u s 

dt 8t 0y S 

For a right-hand coordinate system and constant angular velocity 
S 1 = O* y k 


( 52 ) 


(53) 


Therefore, 

+ € ijk R* u k 

dt at 


(54) 


which is equivalent to the classical result obtained using vector analysis. 

Formulation of the equations of change in a completely arbitrary coordinate system 
is particularly useful in problems which involve moving boundaries, discontinuities, or 
free surfaces. By locating a parametric line along the moving boundary, boundary con- 
ditions become trivial and the computational algorithm is simplified. As an example, 
consider the flow of blood through a single chamber of the heart. This system is typified 



Figure 2 .- Typical open system with moving boundaries. 

in figure 2. This flow problem could be formulated in an Eulerian sense by covering 
the region of interest by an appropriate time invariant coordinate system. Figure 3 
shows a typical approach. However, for fixed grid spacing, numerous problems are 
associated with boundary conditions. These problems are particularly complex when 
the boundary is moving. The problem of boundary conditions could be eliminated by 
adapting a Lagrangian formulation. However, numerous computation problems arise 
in the presence of vorticity which tends to twist the parametric lines into complex over- 
lapping geometries. 
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As an alternate approach, consider a time varying coordinate system for which the 
boundary is along a parametric line. This condition is generated by requiring that on the 
boundary 


n 1 S i = n* 


(55) 


where n 1 is the outward boundary normal and Uj is the continuum velocity at the 
boundary. 

In addition to locating a parametric line along the boundary, the coordinate velocity 
field interior to the boundary must also be constructed. One simple method of defining 
this geometry is through scaling of the boundary surface. Under this scheme, a typical 
representation of the arbitrary coordinate system is shown in figure 4. In this example, 



Figure 4.- Representation of arbitrary coordinate system. 

orthogonality between the chamber normal and scaled chamber tangential lines is 
assumed. 
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An alternate approach for specifying the internal geometry is to generate an 
ir rotational coordinate velocity field. Let the coordinate velocity be defined as the 
gradient of a scalar 

S 1 = g lj *|j (56) 

where <f> is the coordinate velocity potential. The divergence of this vector field is 
given by 



$ 


lij 


The mean divergence of the vector field over the region is given by 


(57) 


D= ^^ n i SldA (58) 

Equation (58) is evaluated in units consistent with the moving coordinate system. Assume 
that the divergence is constant throughout the region; therefore 


or 

gl] *|ij =D (60) 

Equation (60) can be solved for the coordinate velocity potential and hence (using eq. (56)) 
for the coordinate velocity field. The methods described above for generating the coor- 
dinate velocity field are illustrated in appendix B, 

As an additional example, consider the problem of atmospheric flow over complex 
terrain. In order to simplify the boundary conditions, a coordinate system can be defined 
for which the terrain is a parametric surface. For problems that are strongly condi- 
tioned by the time -dependent height of the inversion layer, it may be convenient to define 
this layer as an additional parametric surface of the coordinate system. As a result, the 
coordinate system has time-dependent metrical properties. Scaling between the paramet- 
ric surfaces is arbitrary. A convenient choice is exponential scaling in order that the 
coordinate stretching will approximately parallel the atmospheric density variation. This 
coordinate system is illustrated in figure 5. 
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Time-dependent inversion height 



'<xr 


Terrain 

Figure 5-- A representation of scaled parametric surfaces. 


CONCLUDING REMARKS 

An expression for the co- moving derivative of a tensor has been derived. The 
expression derived is valid for coordinate systems having time-dependent metric ele- 
ments. Time dependency of the metric elements was related to a variable describing 
the coordinate velocity. The resulting theory represents a generalized viewpoint (of 
which the Eulerian and Lagrangian viewpoints are limiting cases) of the field equations of 
continuum mechanics and is useful in problems containing moving boundaries or 
discontinuities. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 31, 1974. 
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APPENDIX A 


A DESCRIPTION OF COVARIANT DIFFERENTIATION 
IN TERMS OF TRANSFORMATION LAWS 

The purpose of this appendix is to show how the covariant derivative can be 
expressed in terms of the elements of the tensor transformation law and their derivatives. 
This development is similar to a related discussion in reference 4. 

Consider a first-order covariant tensor Aj. Let A^x) and A^y) be components 
of the tensor A. in an inertial rectangular Cartesian coordinate system xj and an 
arbitrary coordinate system y i? respectively. Assuming that A^y) has continuous 
first partial derivatives, the covariant derivative is defined as 


The equivalent definition of the covariant derivative of the associated contravariant first- 
order tensor is 


The Christoffel symbol of the second kind can be written as 



where the Christoffel symbol of the first kind is given by 


(A2) 


(A3) 



An element of the metric can be expressed in terms of the elements of the covariant 
transformation law and written as 
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APPENDIX A - Continued 


= T ik T jk 


where 


T ik = 


ax* 

Sy 1 


Assuming continuous first partial derivatives 


0yi ( Tqk ) dy^ ik 


the Christoffel symbol of the first kind becomes 


r i ST. 

[ lS ’ r ] = T rq-^ 


(A5) 


(A6) 


(A7) 


(A8) 


Consider the matrix T formed by the elements TV. Let X- be an element of the 
inverse of this matrix. Then, 


A. T — I 
kr rq 


Note also that, by definition, 


g Pr g- = 6 P 


sr 


(A9) 


(A10) 


By substituting equation (A8) into equation (A3), the Christoffel symbol of the second kind 
can be written as 



(All) 


which, by substituting equations (A9), (A5), and (A10), reduces to 



(A12) 
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APPENDIX A - Concluded 


As a consequence of equation (A9), the Christoffel symbol of the second kind can also be 
written as 



(A13) 


Therefore, for a covariant first-order tensor, the covariant derivative becomes 



9 A i (y) 

8y s 


ax, 


T. 


qp 


iq 


ay £ 


A p (y) 


(A 14) 


An equivalent expression for a contravariant tensor can be deduced from equation (A2). 
The result can be written 



APPENDIX B 


THE CO-MOVING DERIVATIVE OF VELOCITY IN A HARMONICALLY 
PULSATING TWO-DIMENSIONAL POLAR COORDINATE SYSTEM 


Consider the two-dimensional region bounded by a harmonically pulsating circle. 
Over this region, define an orthogonal coordinate system in which the bounding circle is 
a parametric line and the family of similar parametric lines are generated by linear 
scaling. This geometry is shown in the following sketch: 



Assume that the bounding circle has the following harmonic behavior: 


R = R^l + £ sin(o>t - a)J (Bl) 

where the mean radius R, the amplitude of the oscillation 4, the frequency u>, and the 
phase angle a are constants. 

In this simple example, it is possible to write down a coordinate transformation and 
determine the metric elements directly. The coordinate transformations can be written 


xl = 


x 


2 = 


R cos y^ 
(• R sin y^ 

y 1 

V b / J 



(B2) 


where R is a function of time and where y^ is the y* coordinate at the boundary 
(a constant). By definition, the metric elements are 
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APPENDIX B - Continued 


g 12" g 21 _0 1 



j 


Note that these elements are time dependent; that is, 






(B3) 


(B4) 


For a two-dimensional time -dependent polar coordinate system, components of the 
coordinate velocity are 

S 1 » /ziW 
\R 8t \ 

W f (B5) 

S 2 = 0 


In general, the co- moving derivative of the velocity u 1 is 


K) = ^ +v ^ u i +s i u k 

\dt ) »t |<* |k 


(B6) 
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APPENDIX B - Continued 


For a two-dimensional polar coordinate system, all Christoffel symbols vanish except the 
following: 



Therefore, the components of the co- moving derivative become 


and 


duA = y l + v 2/V _ 1 u 2\ 

\U/l 8t 8yl ^y 2 / R 8 1 



8t V yV 


+ v 2 + 

9y 2 y 2 


(B8) 


(B9) 


It is interesting to note that the harmonically pulsating polar system is also gen- 
erated from the assumptions of an irrotational coordinate velocity field and constant 
divergence over the region. The velocity of the boundary is equal to the time derivative 
of the radius transformed to the y s coordinate system. The components are 


.1 


S 1 = 


b 9R 


R 9t 


(BIO) 


S 2 = 0 




Therefore, the mean divergence for the bounded region is given by 


D = 


9 1 ™ 

1 R 9t 


(Bll) 
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APPENDIX B - Continued 


As a result, 


ij - 0 1 9R 

S |lJ R 9t 


(B12) 


The initial values of the metric elements are defined by the initial coordinate geom- 
etry. Assume that, initially, the parametric lines are concentric circles. Therefore, at 
t = 0, 


V 

Sl2 = S 2 l = 0 



^22 



(B13) 


Expanding equation (B12) gives 
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which has the solution 



(y 1 ) 5 
R atV / 


+ Constant 


(B14) 


(B15) 


The corresponding velocity field is 


q1 „11 

S = g 


$1 


which becomes 


S 


1 


v i iM 

3 R at 


(B16) 


(B17) 
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APPENDIX B - Concluded 


Also, 


q 2 22 . 

S =g *| 2 


which becomes 


S 2 = 0 


Note that this velocity field is equivalent to the one originally specified. 
Time variation of the metric elements is given by 


!5u = Si +S) 

at *|j )|t 


which gives 


3g ll = 2 R 8R 
3t / 1 \2 3t 

y b 


8 !p = 2 R l4-\ is 


a «12 S «12 

3t 3t 


= 0 


(B18) 


(B19) 


(B20) 


(B21) 


# 

Since these derivatives are equivalent to the derivatives obtained in equations (B4), the 
coordinate velocity field is identical to the one specified initially. 
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